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Abstract 

We study the quahtative features of the QCD phase diagram in the context 
of the hnear quark-meson model with two flavours, using the exact renormalization 
group. We identify the universality classes of the second-order phase transitions and 
calculate critical exponents. In the absence of explicit chiral-symmetry breaking 
through the current quark masses, we discuss in detail the tricritical point and 
demonstrate that it is linked to the Gaussian fixed point. In its vicinity we study 
the universal crossover between the Gaussian and 0(4) fixed points, and the weak 
first-order phase transitions. In the presence of explicit chiral-symmetry breaking, 
we study in detail the critical endpoint of the line of first-order phase transitions. We 
demonstrate the decoupling of pion fiuctuations and identify the Ising universality 
class as the relevant one for the second-order phase transition. 



1 Introduction 



The phase diagram of QCD at non-zero temperature and baryon density has been the 
subject of many studies during the last years. (For a review see ref. [1].) Parts of this 
phase diagram may be observable through the heavy-ion collision experiments at RHIC 
and LHC In particular, the conditions for the transition from the hadronic phase to 
the quark gluon plasma at high temperature and relatively low baryon density could be 
realized at these experiments. It is, therefore, desirable to connect the various regions of 
the phase diagram with clear experimental signatures. The most prominent feature of the 
phase diagram in this respect is the critical endpoint that marks the end of the line of 
first-order phase transitions. Its exact location and the size of the critical region around it 
determine its relevance for the experiments. There exist several analytical studies of the 
phase diagram near the critical endpoint. They employ effective theories of QCD [2]-[10], 
as analytical calculations in the context of the exact theory are very difficult. Lattice 
simulations are a source of precise information for full QCD. Recent studies [11, 12, 13] 
have overcome the difficulties associated with simulating systems with non-zero chemical 
potential, even though the continuum limit or realistic quark masses have not been reached 
yet. 

The investigation of large regions of the phase diagram, so as to capture the interplay 
between critical and non-critical behaviour around the endpoint, seems more accessible 
to analytical methods. The description of the critical region makes the use of the renor- 
malization group necessary. The study of first-order phase transitions is easier within the 
Wilsonian (or exact) formulation [14]. For these reasons, we would hke to derive the QCD 
phase diagram employing an analytical approach that incorporates the exact renormal- 
ization group. A lot of work has been done on this method in recent years [15]. We shall 
use the effective average action [16], that results by integrating out (quantum or thermal) 
fluctuations above a certain cutoff k. Its dependence on k is given by an exact evolution 
equation [17, 18], while it becomes the effective action for A; = 0. It can by used in order 
to desribe very efficiently the universal and non-universal aspects of second-order phase 
transitions, as well as strong or weak first-order phase transitions [19]. 

We shall not study the fuU QCD for several reasons: the various formulations of the 
exact renormalization group for gauge theories [20]-[23] have not been complemented yet 
with efficient approximations for practical calculations. Moreover, a treatment of the 
exact evolution equation that could capture the effects of the confining region of QCD 
has not been developed. 

We shall work within an effective theory of QCD, the linear quark-meson model, that 
uses as fundamental degrees of freedom mesons coupled to quarks. If only the two lightest 
quarks are taken into account the model is similar to the cr-model of Cell-Mann and Levy 
[24], with the nucleons replaced by quarks. The Langrangian involves four scalar fields 
(the (T-field and the three pions tTj) arranged in a 2 x 2 matrix $, as well as the u and 
d quarks in a doublet. It is invariant under a global SU{2)l x SU{2)r symmetry acting 
on $ and the quark doublet. If the cr-field develops an expectation value the symmetry 
is broken down to SU{2)l^r. The model is very useful for studying the spontaneous 
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breaking of chiral symmetry, using $ as the order parameter. The exphcit breaking of 
this symmetry through the quark masses can also be incorporated. It corresponds to the 
interaction with an external source through a term —ja in the Lagrangian. 

The linear quark-meson model has been studied for a general number of flavours and 
its phenomenological viability has been established [25]. Its formulation for non-zero 
temperature and quark chemical potential has also been developed [26, 27]. In particular, 
it has been found that, for two massless quark flavours, the model predicts a second-order 
phase transition for increasing temperature and zero chemical potential, and a flrst-order 
one for increasing chemical potential and zero temperature. It is likely, therefore, that 
the expected phase diagram of QCD for small chemical potential can be reproduced 
correctly within this effective theory. We point out, however, that the study of the colour 
superconducting or colour-flavour locked phases [1] requires the introduction of degrees 
of freedom that correspond to diquark condensates. 

In this work we shall study the phase diagram of the quark-meson model. We consider 
only the case of two flavours as the basic structure appears already within this approxi- 
mation. The generalization to three flavours is straightforward, even though the degree 
of complexity increases significantly. We do not consider realistic values for the various 
parameters of the Lagrangian. The reason is that the treatment of the exact evolution 
equation becomes complicated for the realistic values and little analytical progress can be 
made. Instead, one has to rely heavily on numerics that obscurs the basic features of the 
problem. The consideration of parameter values consistent with the phenomenology is an 
obvious subject for future studies. 

Despite the above limitations, our work includes several important elements: The 
infrared region is treated correctly through the renormalization group. The fixed point 
structure can be obtained in detail. The universal behaviour near the critical points can 
be distinguished from the the non-universal aspects. The various universality classes can 
be identified and the crossover behaviour from one to the other can be studied. 

In section 2 we discuss the phase diagram in the context of mean field theory in order 
to determine the main features we expect to appear. In section 3 we present the main 
ingredients of the model we are considering. In section 4 we derive the evolution equation 
that describes the dependence of the potential on a scale k that acts as an infrared cutoff. 
In section 5 we cast this equation in a form that has no explicit dependence on k and 
facilitates the identification of fixed points. In section 6 we derive the initial conditions 
for the solution of the evolution equation at energy scales below the temperature, where 
the evolution becomes effectively three-dimensional. In section 7 we use the formalism in 
order to discuss the phase diagram in the absence of explicit chiral symmetry breaking 
(zero current quark masses). In section 8 we discuss the phase diagram in the presence 
of explicit chiral symmetry breaking. In section 9 we summarize our findings and outline 
the directions of future research. 
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2 Mean field theory 



In order to obtain some intuition on the form of the phase diagram, it is instructive 
to discuss a simple model in the mean-field approximation. The Lagrangian contains 
standard kinetic terms for the fields a, iTi {i = 1, 2, 3), and the potential 



with p — ((T^ + 7ri7r*)/2. The couphngs m^, A, u are taken to be functions of two external 
parameters, which we denote by T and /i without having to be more specific at this stage. 
We assume that u is always positive, so that the potential is bounded from below. 

Let us consider first the case A > 0. For > the minimum of the potential is 
located at cr = tTj = and the system is invariant under an 0(4) symmetry. For < 
the minimum moves away from the origin. Without loss of generality we take it along 
the a axis. The symmetry is broken down to 0(3). The vr's, which play the role of the 
Goldstone fields, become massless at the minimum. If m'^{T,fi) has a zero at a certain 
value T = T^r, while X{Tcr,n) > 0, the system undergoes a second-order phase transition 
at this point. The minimum of the potential behaves as ctq ~ |T — T^r]^/'^ slightly below 
the critical temperature. The critical exponent P takes the mean field value P — 1/2. 

Let us consider now the case A < 0. It is easy to check that, as increases from 
negative to positive values (through an increase of T for example), the system undergoes 
a first-order phase transition. For |A| approaching zero the phase transition becomes 
progressively weaker: The discontinuity in the order parameter (the value of a at the 
minimum) approaches zero. For A = the phase transition becomes second order. The 
minimum of the potential behaves as ctq ~ \T — |^/^. The critical exponent (3 takes the 
value /? = 1/4 for this particular point. 

Let us assume now for simplicity that A is a decreasing function of p only, and has a 
zero at /I — /i^. If we consider the phase transitions for increasing T and fixed /i we find 
a line of second-order phase transitions for /i < /i^, and a line of first-order transitions for 
II > fi^. The two lines meet at the special point (T*, /x*), where = T^^ for /i — /i^. This 
point is characterized tricritical point. 

If a source term —ja is added to the potential of eq. (1) the 0(4) symmetry is explicitly 
broken. The phase diagram is modified significantly even for small j. The second-order 
phase transitions, observed for A > 0, disappear. The reason is that the minimum of the 
potential is always at a value a ^ that moves close to zero for increasing T. For small 
j the mass term at the minimum approaches zero at a value of T near what we defined 
as Tcr for A = 0. However, no genuine phase transition appears. Instead we observe an 
analytical crossover. 

The fine of first-order transitions persists for j 7^ 0. However, for increasing T the 
minimum jumps discontinuously between two non-zero values of a and never becomes 
zero. The line ends at a new special point, whose nature can be examined by considering 
the (T-derivative of the potential U{p) of eq. (1). For a first-order phase transition to 
occur, dU/da must become equal to j for three non-zero values of a. This requires 



U{p) =m^p+-Xp' + -up^, 
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d^U/da"^ = at two values of a. The endpoint corresponds to the situation that all 
these values merge to one point. It can be checked that this requires A < 0. At the 
minimum cr* of the potential U (cr; j) — U {a, tTj = 0) —ja at the critical endpoint, we have: 
dU/da = d'^U/da'^ = d^U /da^ = 0. Therefore, at the endpoint wc expect a second-order 
phase transition for the deviation of a from o"*. As A 7^ 0, the critical exponent /? takes 
the value (3 = 1/2. The tt's are massive, because = d'^U{p)/dTif = dU{p)/dp = j/a^ 
at the critical endpoint. 

The purpose of the following is to recontruct the above picture in the context of 
the linear quark-meson model. We shall take into account the effect of fluctuations of the 
fields, so that the potential will have a more complicated form than the simple assumption 
of eq. (1). Moreover, the nature of the fixed points will differ from the predictions of 
mean field theory. 

3 The model 

We consider an effective action of the form [25, 26, 27] 



Tk = J^^"^ dx' J d''xUilj''{Yd^ + ii^')ilJa + hilj'' 



with p = Tr (^$1^$^. The field $ includes the a field and the three pion fields. It is 
parametrized as 

$ = 1 ((J + iTT • f) (3) 

where r; (/ = 1,2,3) denote the Pauli matrices. The fermionic field includes two 
fiavours (a = 1,2) corresponding to the up and down quarks. The non-zero temperature 
effects are taken into account by restricting the time integration within the finite interval 
[0, 1/T] in Euclidean space. The scalar (fermionic) fields obey periodic (anti-periodic) 
boundary conditions. We are also considering a non-zero chemical potential, associated 
to the quark number density ijjlijj"' that is proportional to the baryon number. 

In the effective action we have omitted the wavefunction renormalization for both the 
scalar and fermionic fields. We have also neglected the scale dependence of the Yukawa 
coupling h. We take into account only the scale dependence of the potential Uk, which 
we assume to be a general function of p. Our approximation is consistent in two cases: 
a) If the Yukawa coupling is small at an initial large momentum scale A, its subsequent 
evolution for A; < A is expected to be slow. The reason is that the leading contribution 
to dh? /dink is ~ h?' [27]. The same holds true for the wavefunction renormahzations 
[27]. The scalar field wavefunction renormalization receives an additional contribution 
through scalar field contributions. Typically this is small in three and four dimensions, 
as is obvious from the smallness of the anomalous dimension in the pure scalar theory. 
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2) If partial infrared fixed points exist in tlie evolution equations, all the above quantities 
take almost constant values. The wavefunction renormalizations can be set equal to 1 
through the rescahng of the fields, while the Yukawa coupling becomes scale independent. 
This scenario is realized in the low energy regime of the quark-meson model in four 
dimensions, for realistic values of the quark and meson masses. 

The last term in eq. (2) accounts for the explicit breaking of the chiral symmetry 
through the current quark masses. One expects j ~ M = diag (mu,md). We assume 
equal masses m for the two light quarks so that ~ mS°'b- The determination of the 
proportionality constant requires the embedding of the theory described by the action (2) 
in a more fundamental framework. For example, in refs. [27, 28, 29] it was calculated 
in terms of the energy scale at which $ can be introduced as a composite field in 
a theory of quarks (such as the Nambu-Jona-Lasinio model [30]). One obtains j — 
2Z^^f;^ml^M/hk^, with Z^^fe* the fermion wavefunction renormalization, and the 
positive mass term in the potential of $ at the scale k^. This mass term becomes negative 
at lower energy scales because of fcrmionic quantum fiuctuations, and the chiral symmetry 
is spontaneously broken. In this work wc set Z.ip^k = 1 for all k and also neglect the running 
of /ifc. For this reason it is not possible to express j in terms of fundamental quantities. 
However, j can be fixed through the low energy dynamics. The pion mass at the vacuum 
in the presence of explicit chiral symmetry breaking is proportional to the source term 
and can be used as input for its determination. 



4 The evolution equation 

The dependence of the effective action on an infrared cutoff scale k is given by an exact 
fiow equation [17], which for fermionic fields ip (quarks) and bosonic fields $ (mesons) 
reads 



f^m,^ = iiv, (rf 1^,*] + - TV, (rf 1^,*] + . 

(4) 

(2) 

Here is the matrix of second functional derivatives of Fj. with respect to both fermionic 
and bosonic field components. The first trace in the r.h.s. of (4) runs only over the 
bosonic degrees of freedom. It implies a momentum integration and a summation over 
fiavor indices. The second trace runs over the fermionic degrees of freedom and contains 
in addition a summation over Dirac and color indices. The terms RkB, RkF are infrared 
cutoffs that have been introduced in order to eliminate contributions from fiuctuations 
with momenta ^ k"^ in the momentum integrations. 

The evolution equation for the scale dependent potential Uk has been computed in 
ref. [26] starting from eq. (4). In the approximation that we neglect the running of the 
Yukawa coupling it takes the simplified form: 

^C/fe(p; T, //) = ^UkB{p] T, /i) + ^UkFip; T, //) , (5) 
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where 



d 

—UkB{p;T,ii) 



2 ^ 



(27r)3 dk 



+ 



Pk{q' + 4nW^) + U',{p;T,i^) 

(6) 



and 



d_ 
dk 



d_ 
dk 



(7) 



1 ^ 



-AN^T 



-oo 
oo 

I 



(27r)= 



In 1 + exp 



g2 + A;2 + /i2p/2 - P 



.(9) 



Here p = Tr = (cr^ + tt^) /2 and Nc — S denotes the number of colours. 

The two terms in the r.h.s. of eq. (5) correspond to the contributions from scalar 
and fermionic fluctuations. The scalar term, given by eq. (6), includes the contributions 
from the radial mode (the a field) and the three Godstone modes (the pions). The 
parametrization in terms of p makes apparent the 0(4) symmetry characterizing the 
mesonic sector in the two-flavour case. The meson mass terms are U^. + 2pU'^, U'^, where 
primes denote derivatives of the potential with respect to p. If the chiral symmetry is 
broken by an expectation value for the a field, the mass terms become = d^Uk/da'^, 
ml- = {1/ a)dUk/ da , as expected. 

The scale dependent propagators also contain the momentum dependent part Pk = 
+ RkB{q^)- The term 



RkB{q^) 



9 -„2//,-2 

1 _ p-«Vfe^ 



(10) 



acts as an effective infrared cutoff in the loop integrations that determine the effective 
potential. The scale dependence of Uks originates in this cutoff, as is apparent from the 
/c-derivative of RkB in eq. (6). The periodic conditions imposed on the scalar fields for 
non-zero temperature result in the replacement of the integration over g° by a sum over 
the discrete Matsubara frequences. The chemical potential p affects U^b only through 
the contribution of the fermionic part U^f to the mass terms. 

The fermionic contribution Ukp has a simple interpretation within our approximate 
treatment. It consists of three parts: The term '^VkF^p) originates in the quantum fiuctu- 
ations of the fermionic fields, for which a mass-hke infrared cutoff A;^ has been introduced. 
This type of cutoff does not provide automatic ultraviolet regularization. This is achieved 
through the modification of the cutoff by an additional ^^-function [26]. The number of 
degrees of freedom is 2(flavour) x2(spin) x2(particle-antiparticle) x A^c (colour). The term 
^ykF{p',T,iJ,) is the free energy of a non-interacting fermionic gas of temperature T and 
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chemical potential /j, in the presence of the mass-like cutoff k^. The ultraviolet regular- 
ization has been omitted as we assume A » T. The number of degrees of freedom is 
2 (flavour) X 2 (spin) xA^c (colour). The term ^Vkpip^T, —jj) is the free energy of the anti- 
particlcs that have opposite chemical potential. We emphasize that the evolution of Ukp 
is derived starting from the exact flow equation for the effective action of the system. The 
truncation of eq. (2) and the additional assumptions of constant wavefunction renormal- 
izations and Yukawa coupling lead to the intuitive form of eqs. (7)-(9). We also point out 
that the use of different cutoffs for the scalar and fermionic fields does not affect physical 
quantities, as the cutoffs are removed in the limit A; ^ 0. 



5 Scaling form of the evolution equation 

Rather than solving eq. (5), it is more convenient to solve the evolution equation for Vy.. 
This can be written in the form 

^Ul{p;T,t,) = ^4A;(3C/^ + 2pC/ni^t((C/^ + 2pC/DMT)+3^4^C/^i^f(c/^A^T) 



+ 



'VkF{p) + 'VkF{p;T,p) + 'VkF{p;T,-p) , (11) 



dkdp 

where, in general dimensions d, 



vt - s'^+V'^/^r (A . (12) 



The "threshold" functions L\ are particular cases of 

L^(w;T) = -2nA;2"-'^+i7r-i+^r (-\ T^, [ d'^-^q ^{Pk + w kY^^+'\ (13) 

\2J ^ J ok 

where is a function of (f + An^Tr'^T'^, as we discussed earlier. 

The functions L'^{w;T) have been discussed extensively in refs. [31, 32, 18, 19]. Their 
basic properties can be estabhshed analytically. For T <^k the summation over discrete 
values of m in the expression (13) is equal to the integration over a continuous range of 
(f, up to exponentially small corrections. Therefore 

Liiw- T) = Li{w, 0) = Li{w) for T « k. (14) 

In the opposite limit T ^ k the summation over m is dominated by the m = contribu- 
tion. This results in the simple expression 

Li(w;T) = ^-^^Lt\w) for T » k, (15) 

with Vd defined in eq. (12). The two regions of T/k in which L'^{w;T) is given by 
the equations (14), (15) are connected by a small interval < k < T/9i with a 
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more complicated dependence on w and T. The transitions a.t k = T/9i and k = T/62 
are sharp for an exponential cutoff, such as the one given by eq. (10). The functions 
L^(u') fall off for large values of w, following a power law. As a result they introduce a 
threshold behaviour for the contributions of massive modes to the evolution equations. 
These include LJ^ functions with the mass eigenvalues divided by k"^ as their arguments. 
When the running squared mass of a massive mode becomes much larger than the scale 
/c^, the mode decouples the respective contribution vanishes. We evaluate the expressions 
for T) numerically and use numerical fits for the solution of the evolution equations. 

The fcrmionic contribution is expected to display a similar behaviour for decreasing 
k. It can be checked that the temperature-dependent term of eq. (9) evolves slowly 
for A; ^ T and obtains asymptotically a constant form. However, the k dependence 
disappears only as power law, instead of exponentially fast as for L'^{w,T) above. The 
quantum contribution of eq. (8) has a strong dependence on the ultraviolet cutoff A that 
regulates the momentum integration. As a result, the quantity {^Vkp) / dkdp displays 
a logarithmic dependence on A. The above shortcomings are connected to the use of 
a mass-like cutoff for the fermionic sector, instead of an exponential cutoff as for the 
bosons. A mass like cutoff is the simplest choice that preserves the Lorentz structure of 
the kinetic term of free fermions for non-zero chemical potential, so as to make practical 
calculations feasible [26]. A more ingenuous, and more complicated, form is required in 
order to preserve the correct Lorentz structure and guarantee exponential decoupling of 
the temperature effects for fermions at low energy scales. 

In this work we are not interested in the details of fermion decoupling. For this 
reason we follow a different approach. We assume that at a sufficiently low scale k ^ T 
the fermionic contributions to the evolution equation decouple. This is the expected 
behaviour, as the lowest Matsubara frequency is ttT for the fermions and no zero mode 
exists. The result of the integration of the fermionic modes is a contribution equal to the 
perturbative one-loop term for the potential. It is clear that this assumption is consistent 
with the form of the terms of eqs. (8), (9) for k — 0. It is expected to be correct for small 
Yukawa coupling h. 

We concentrate on the behaviour of the evolution equation (11) for k ^ T, in which, 
as we explained above, two important simplifications occur: 

a) The "threshold" functions are approximated very well by Lf{w,T) = {4:T/k)L{{w). 

b) The fermions have decoupled because of their thermal masses. 
We define the quantities 



(16) 



T 

P_ 



(17) 



and their dimensionless versions 



Uk{p) 



fc3 



k^T 



(18) 



P 



El = JL 

k kT 



(19) 
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The evolution equation (11) can be rewritten in the scahng form [18] 



-u'M = -2< + ~pul + (3< + 2pO Ll {u'k + 2~pul) + 3t;3 ul «) , (20) 



with t = ln(/c/M) and the primes denoting derivatives with respect to p. The energy scale 
M can be defined arbitrarily. The characteristic property of the above equation is that 
the scale k does not appear explicitly in it. This permits the search for fixed points, i.e. 
scale independent solutions. They are the solutions with du'^/dt = 0. 

The evolution equation (20) has a form typical of a three-dimensional theory (dimen- 
sional reduction). This is expected, as only the fields with zero Matsubara frequences 
(the scalar zero modes) are sufficiently light in order to contribute to the evolution. 

6 Initial conditions for the evolution equation 

The initial conditions for the solution of eq. (20) must be given at the scale = T/92- 
Let us assume that at the large momentum scale A (where the effects associated with the 
temperature and chemical potential are negligible) the potential has the simple form 



with poA < 0. For small A wc can integrate the evolution equation (11) between the scales 
A and kr = T/62, approximating the arguments of the threshold functions with zero and 
setting U'^ — U'^' — 0. These approximations account for the leading pertubative result 
and give [32, 33] 



with ^VoF, "^Vqf given by cqs. (8), (9) with k = 0. We have assumed that the fermionic 
contributions have already decoupled at kx = T/62 and are given by eq. (9) with A: = 0. 
The term ~ results from the complicated form of L\{0] T) in the /c-interval [T/ ^2, T / 9i\ 



c/a(p) = c/a(p;0,o) = -(p-poa), 



(21) 




(22) 



[33]. 



The T — fj, — fermionic contribution (8) can be evaluated exphcitly 




1 + 



2A2 



2A^ 



(23) 



We can cast eq. (22) in the form 




02 




+ If{p;T,p) 



(24) 
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with 



Pgr — PoA 



lF{p;T,fj,) 



167r2 



A2 + 



167r2 



pin 



A 87r2 

1 -\- ^ 



A' (25) 
Vof(p; t, /.) + -//). (26) 



The T — II — potential displays the spontaneous chiral symmetry breaking induced 
by the fcrmionic fluctuations. For sufiiciently large Yukawa coupling h, even if poA is 
negative (corresponding to a positive mass term), the renormalized value Pqr can become 
positive through the fermionic contribution in eq. (25). 

At non-zero temperature, a consistency check for the initial condition (24) can be 
performed by integrating eq. (20) from k — T/92 to k — 0. Again we assume a small Xr, 
set U'^ — Xr, U'H — and take the arguments of the "threshold" functions equal to zero. 
In this way we obtain 



f/^(p;T,p)~t/' (p;r,p) + 



47r2 62 



XT^ = X 



P- PoR + Q 



rp2 

24 



+ /p(p;T,p), (27) 



where we have used ^f(O) = —\pn. This expression reproduces correctly the leading 
perturbative result [34] for the critical temperature. However, the assumption that the 
couplings remain small is not valid near the phase transiton. Close to the critical tem- 
perature the fixed-point structure of the theory becomes important and the arguments of 
the "threshold" functions cannot be neglected. 

Another ingredient that we have neglected in our considerations so far is the explicit 
chiral symmetry breaking term ~ cr, arising through the current quark masses. A term 
linear in a does not appear expicitly in the evolution equations, as the flow equation (4) 
involves second functional derivatives of the action with respect to the flelds. For this 
reason the linear term can be added directly to the potential at any scale. 



7 The phase diagram for j = 

Let us summarize the main points of the framework we have set up. Our theory is defined 
through the potential of eq. (21) at a scale A » T. In a reahstic theory we can identify A 
with the energy scale /cj, at which the mesonic fields a, tTj can be introduced as composite 
fields in a theory of quarks. The renormalized potential for A; = and zero temperature 
and chemical potential is given by eqs. (24)-(26) with T = p = 0. These approximate 
expressions are correct if: 

a) h is small enough for its logarithmic running to be neglected. This requires corrections 
~ /i^/(47r^) to be negligible relative to Ir? . 

b) A is small enough for terms ~ A^ to be neglected. Notice, however, that corrections 
~ I (47r^) arc not assumed to be small relative to A. On the contrary, it is for the range 
h'^ I (47r^) ~ A that first-order phase transitions appear. 
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The explicit chiral symmetry breaking can be taken into account by adding a term —ja 
to the potentiaL The parameters Pqr, A, j can be determined through the pion decay 
constant and the masses rrifj, m-^ of the a and pion fields. The Yukawa coupling h can 
be determined through the constituent quark mass My. 

At the scale /ct = T/62 the potential is given by eq. (22), with ^'Vop given by eq. (9) 
for k = 0. The evolution of Uk for k < kx is given by eq. (20). The scalar modes with 
non-zero Matsubara frequencies, as well as the fermions, have completely decoupled. Only 
the fiuctuations of the scalar zero modes contribute to the evolution of the potential. For 
this reason the theory is effectively three-dimensional at these low scales. For practical 
calculations we use ^^2 = 1- 

The fermionic contribution Ji?(p;T, /x) to the initial condition for the potential, given 
by eq. (26), can be evaluated in certain limiting cases for small mass term h'^p/2. For 
II — the standard high temperature expansion [34] gives 

I,(p; T, 0) = AT.^T^ + iV.^pln + 0{h% (28) 

with a ~ 8.5. For T = the Fermi-Dirac distribution becomes a step-function and we 
obtain 

I4p;0,p) = N^^y + iV.^pln + Oih% (29) 

with (3 ~ 10.9. 

The terms ~ h^p in the r.h.s. of the above equations can change the sign of the 
tree-level term Xp in the initial condition (24). The phase diagram expected for QCD 
is obtained if the total term ~ p is positive near the phase transition for T ^ 0, p — 
and negative for T = 0, // 7^ 0. In our model we choose couplings such that this 
condition is satisfied. It has been checked explicitly [26, 27] that the realistic values of 
the parameters, consistent with low-energy QCD phenomenology, reproduce the correct 
order for the transitions for T 7^ 0, /i = and T = 0, ;U 7^ 0. 

We define dimensionful quantities in terms of the parameter poi? defined in eq. (25). 
Our analysis of the phase diagram will be carried out for a model with A = 0.1, h'^ — 2.3, 
A2 = 1.8poi?- 

7.1 The second-order phase transitions 

We start by considering a system with zero chemical potential. In fig. 1 we display the 
solution of eq. (20), for a theory with p = 0, T / ^ypoR — 0.407. Wc start the evolution 
at a scale kr = T and move towards the infrared. The initial condition is given by eqs. 
(24), (26), and corresponds to a system in the broken phase. We observe that initially 
u'f^{p) is small, because of the smallness of the value of A we chose. Subsequently the 
potential evolves (dashed lines), while its minimum po remains non-zero. At a later stage 
the solution approaches an almost constant form (solid lines). This is the fixed point of 
the 0(4) imiversality class. At the last stages of the evolution (dotted lines) the minimum 
of the potential moves to the origin and u'^^iQ) becomes positive. 
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Figure 1: The evolution of the rescaled potential for fj. = 0, T/ ^/poR — 0.407. 

For A; — > 0, the rescaled mass term ^'^.(O) diverges so that (i[4(0; T)/dp = k'^duk{0)/dp 
approaches a constant value. The system ends up in the symmetric phase. The renormal- 
ized mass term m\{T) = dUo{0;T)/dp is very small in terms of the characteristic mass 

1 /2 

scale of the potential at k = T. The reason is that, during the "time" t = ln(A;/pQ^) 
the solution is close the fixed point, the quantity 'u'^(O) remains constant, so that U'j,{0; T) 
scales ~ k'^. The longer the system spends near the fixed point, the smaller the final 
value of mj^{T). The approach to the fixed point is controlled by one parameter, the 
temperature T. For the evolution displayed in fig. 1 the temperature is slightly larger 
than the critical value T„. that separates the symmetric from the broken phase. For a 
value slightly smaller than the critical, the evolution would be similar to that in fig. 1, 
apart from the last stages. For k ^ the minimum po of Uk{p) would diverge, so that the 
minimum po of Uk{p;T), given by pQ{k,T) = p{k)kT, would reach a constant very small 
value por{T). In this way, the system would end up in the broken phase. Examples of 
this behaviour can be found in the review [19] and references therein. 

The fixed point of the 0(4) universality class is approached by fine-tuning only one 
relevant parameter, the temperature T of the system. The temperature effectively controls 
the mass term of the initial potential. The properties of the critical theory (the one for 
T ~ Tcr) are determined only by the fixed point. The reason is that the system, after 
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Figure 2: The critical exponents v and (3 as the critical temperature is approached. 



spending a significant part of the evolution near the fixed point, loses memory of the 
initial conditions. As a result, all physical quantities at the vacuum can be parametrized 
in terms of the deviation of T from Tcr- The dependence on T — Tcr is not analytic. In 
particular, for the mass term and the expectation value we expect a behaviour that can 
be parametrized as 



Pm{T) 



\T - 71, 



\2v 



1 2/3 



(30) 
(31) 



In fig. 2 we plot the effective values of the critical exponents i/, and /3, extracted through 
the integration of the evolution equation for various values of T. We observe that the 
exponents become constant very close to the critical temperature, with values v ~ 0.82, 
13 ~ 0.41. These are in reasonable agreement with the values determined through high 
precision calculations using alternative methods, or the effective average action taking 
into account the wavefunction renormahzation: v ~ 0.74, 13 ~ 0.38 [19]. The descrepancy 
arises because we have neglected the wavefunction renormahzation, effectively using 77 = 
for the anomalous dimension of the field. The correct value of this quantity is 77 ~ 0.036. 
Notice, however, that the scaling law (3 — r\)v is satisfied by our result. 
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Figure 3: The evolution of the rescaled potential for n/ ^Jpqr — 0.735, T / ^/pm — 0.1. 
7.2 The first-order phase transitions 

We turn now to the region of first-order phase transitions, which are expected to take 
place for large values of the chemical potential. We consider a system with // ~ 0.735, 
T / y/poR — 0.1. In fig. 3 we display the evolution of u'^.{p) starting at A; = T and moving 
towards k = 0. This figure should be compared to fig. 1. The differences arc significant. 
Already at k = T, there are two values of the order parameter p for which m'^ vanishes. 
Moreover, ^'^-(O) is positive. This implies a potential with two minima (one located at the 
origin) and a maximum. During the evolution towards the infrared the system does not 
approach a fixed-point solution. For A; — > 0, u'^{0) increases so that U'f^{0;T, /j.) reaches a 
constant value. The point po where the second minimum of Uk is located diverges also, 
so that the second minimum po of Uk approaches a constant value. 

The behaviour of Uk{p',T,p) is displayed in fig. 4. Wc observe the two minima, and 
the curvature of the potential around them that determines the mass terms, settling down 
at constant values. The minimum at the origin is slightly deeper than the one at non-zero 
p, indicating that the temperature is slightly larger than the critical one. The barrier 
separating the two minima remains scale dependent even after the minima have become 
scale independent to a good approximation. This is a general feature of the solutions 
near first-order phase transitions, which is related to the issue of the convexity of the 
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Figure 4: The evolution of the potential for /x/ ^/poR — 0.735, T/ ^/poR ~ 0.1. 

effective potential [31, 19]. As the scale-dependent potential Uk approaches the effective 
potential for /c — > 0, the barrier should disappear, giving its place to the standard Maxwell 
construction. This is a general property of the solutions of the evolution equation [31, 19], 
which can be established analytically. However, the numerical integration of the evolution 
equation in this regime is difficult because of the presence of a pole in the threshold 
function. For this reason, in fig. 4 we have stopped the evolution at a non-zero value 
of k, such that the minima are stabilized, even though the barrier has not disappeared 
yet. An example of a more elaborate numerical integration that displays the approach to 
convexity for /c ^ is presented in rcf. [19]. The calculation of bubble-nucleation rates 
during a first-order phase transition is performed at non-zero values of k. The physical 
value of the nucleation rate is independent of the choice of k at which the calculation is 
carried out [35] . 

7.3 The tricritical point 

The second-order phase transitions, discussed in subsection 7.1, form a fine on the phase 
diagram, which starts at /i — and continues towards non-zero values of /i. Similarly, the 
first-order phase transitions, discussed in subsection 7.2, form another line that starts at 
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Figure 5: The evolution of the rescaled potential for /// ^JpQR — 0.718, T j ^Jpqr — 0.130. 

T = and continutes towards non-zero values of T. These two lines are expected to meet 
at a point on the (T,//) plane, characterized as the tricritical point. Near it the first-order 
transitions are very weak (the discontinuity in p approaches zero) and eventually become 
second order. 

In fig. 5 we display the evolution of m'^(p) for p/ ^Jpqr = 0.718, T/ ^poR — 0.130. This 
figure is very similar to fig. 1. The potential evolves towards the 0(4) fixed point, where 
it becomes stationary for a while. Eventually it moves away from the fixed point and 
settles down in the symmetric phase. 

The significant difference between the two figures concerns the initial stage of the 
evolution. In fig. 1 this is rather short, as the potential moves quickly away from its 
initial form and approaches the fixed point. In fig. 5 the initial stage is much longer, as 
can be deduced from the accumulation of curves on the horizontal axis. In fact, two fixed 
points are apparent in fig. 5: the one characterizing the 0(4) universality class, and the 
Gaussian fixed point which corresponds to u'^{p) = 0. The system starts very close to 
the Gaussian fixed point. This is unstable, so that eventually the system approaches the 
0(4) fixed point. In order to start near the Gaussian fixed point, two relevant parameters 
must be fine-tuned: the temperature and the chemical potential. They control the mass 
term (the term ~ p) and the quartic coupling (the term ~ p^) of the initial potential. In 
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Figure 6: The critical exponents v and (3 as the critical temperature is approached near 
the tricritical point. 



particular, the initial values of these two terms arc very close to zero, so that the term 
~ p"^ (in an expansion of Ukj,{p) around the origin) is important. 

These findings are in very good agreement with the discussion of section 2, which 
was based on mean field theory. The reason is that the Gaussian fixed point results in 
mean-field behaviour, as the effect of fluctuations is negligible because of the smallness of 
the couplings. Our results indicate that the tricritical point must be connected with the 
Gaussian fixed point of fig. 5. The number of relevant parameters is correctly predicted 
to be 2. 

The effective critical exponents {3 can be computed from eqs. (30), (31), similarly 
to subsection 7.1, by keeping ji fixed and varying T. The results are depicted in fig. 6. 
We observe that for \T — Tcrl/y^poR = exp(— 5) - exp(— 8) they take the values (3 = 0.25, 
V = 0.5. These are the mean-field predictions that we derived in section 2. For this range 
of temperatures the potential never approaches the 0(4) fixed point. It stays close to the 
Gaussian one, before moving directly towards the broken or the symmetric phase. For 
temperatures closer to Tcr, the system feels the attraction of the 0(4) fixed point. For T 
very close to Tcr the potential spends a large part of its evolution near this fixed point. 
As a result it loses memory of its initial form near the Gaussian fixed point. The critical 
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Figure 7: The evolution of the rescaled potential for /// ^JpQR — 0.719, T j 

exponents take the values characteristic of the 0(4) universahty class. 

The curves in fig. 6 between \r — Tcr\/'s/pm — exp(— 9) and \T — Tcr\/'\/pm. — 
exp(— 30) are typical examples of crossover curves. They describe the variation of universal 
quantities, such as the critical exponents, as the relative influence of two fixed points 
changes. 

We can confirm our identification of the tricritical point by studying the first-order 
phase transitions near the Gaussian fixed point. In fig. 7 we display the evolution of 
u'^{p) for a theory with jij ^poR = 0.719, T/ y^po/i' — 0.128. The slight increase of the 
value of the chemical potential relative to fig. 5 modifies drastically the evolution. The 
potential stays very close to the Gaussian fixed point for a long "time" t = \n{k/pl/^), 
but eventually moves away from it. (Notice the difference by a factor 100 of the scale of 
the vertical axis between figs. 3 and 7.) Its form is characteristic of a first-order phase 
transition as we discussed in detail in subsection 7.2. However, the transition is very weak, 
as the discontinuity in the order parameter is orders of magnitude smaller than Pqr. 

The tricritical point is located on the plane {p/y/pm, T/^pofi) in the region between 
the points (0.718, 0.130) and (0.719, 0.128). For the critical values of p and Tcr the 
solution of the evolution equation stays close to the Gaussian fixed point for an infinitely 
long "time" t. As a result the critical behaviour is determined completely by mean-field 
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theory. For example, the critical exponent (3 stays equal to 0.25 arbitrarily close to the 
critical temperature. 



8 The phase diagram for j ^ 

8.1 The analytical crossover and the first-order phase transi- 
tions 

The effect of a source term —ja on a second-order phase transition has been studied 
extensively in the statistical physics and field theory literature. It is well known that the 
phase transition disappears and is replaced by an analytical crossover. No singularities 
appear in the physical quantities. Instead these arc described by smooth functions con- 
necting the two "phases". The reason can be understood by considering the potential for 
a constant value of /i, corresponding to the part of the phase diagram where a second- 
order phase transition is expected for j = 0. The vacuum of the theory is determined by 
the relation dUo/da — j. For j 7^ the vacuum is located away from the origin for all 
values of T. If we take a as the order parameter, no signal of a phase transition exists. 
Moreover, both the cr-field and the pions are massive at the vacuum. This implies that 
no significant fluctuations exist in the deep infrared, as all fields decouple at sufficiently 
low energy scales. Consequently, no singularities are expected in physical quantities. 

For very small j, the influence of the fixed point leads to a universal form of the 
potential in the critical region. It is convenient to define the quantity x — ST/a^^^, with 
ST — T — Tcr and P a critical exponent. The potential can be expressed as 

/W, (32) 



da 



with 6 another critical exponent and dUo/da — j [36]. The form of the function f{x), 
characterized as the universal equation of state, is completely specified by the universality 
class. For small j, the temperature dependence of the a-field and pion masses near the 
critical temperature can be extracted from eq. (32) [37] . 

We shall not determine the equation of state of the 0(4) universahty class here, as it 
has been calculated in ref . [27] . There, it is also demonstrated how the interplay between 
universal and non-universal behaviour can be described efficiently by the effective average 
action. Even though the study of ref. [27] is carried out for /i = 0, it displays all the 
characteristic features associated with theories in the small-/^ region of the phase diagram. 

The first-order phase transitions for large jj, that we discussed in subsection 7.2 are 
not affected significantly by the presence of a source term —ja. The minimum at the 
origin moves to a non-zero value of a, as in the case of second-order phase transitions. 
The linear term also affects the relative depth of the two minima of the potential and 
chages the value of the critical temperature. The most characteristic feature of the phase 
diagram for j 7^ is the presence of a critical endpoint that marks the end of the line 
of first-order phase transitions. At this point a second-order phase transition is expected 
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to take place. In the following subsection we discuss in detail the nature of the endpoint 
and determine the relevant universality class. 



8.2 The critical endpoint 

In order to discuss the critical endpoint we need to adapt our formalism to the situation of 
explicit symmetry breaking through the addition of the source term —ja to the potential. 
In particular, we have to define an appropriate order parameter that can describe the 
second-order phase transition. For this reason we expand the field a around a non-zero 
value (7j, so that a — aj + S with (5 <S cr,-. We define a new potential 



Vki5; T, n) = Ukiaj + S; T, //) - j {aj + 5) 



(33) 



and the dimensionless quantities 

S 

The evolution equation takes the form 



kT 



(34) 
(35) 



/k-'T 



'kT 



+ 5 



~\ 2 



/kJ'T 



/kT 



6 



(36) 

where the primes now denote derivatives with respect to 6. 

The last term arises through the pion fluctuations. This is apparent if one considers the 
argument of the threshold function, which is the pion mass term (l/a){dUk/da) divided 
by k^. The term introduces explicit scale dependence in the equation. We expect the 
contributions that include inverse powers of k to become dominant in the infrared. The 
potential u'j^{S) could also include terms that scale with inverse powers of k. We isolate 
the most singular contribution by parametrizing u'f.{5) as 



Ck 



a/Ft' 



(37) 



In the following we shall show that v'f^{6) can reach a scale-invariant form during the 
evolution, so that Vk{S) <^ j/ \/WT for fixed 5 and A; — > 0. Under this assumption, we can 
write eq. (36) as 



d_ 
di 



d 

dk'' ^ 



-3v:- 



U + Ck)T 3 

12 -^1 



i+ck 



(38) 
(39) 
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These equations have a clear physical interpretation: For sufficiently low k, the pions 
(the Goldstone modes) decouple from the evolution of the cr-field (the radial mode) . The 
reason is that the pions are massive in the presence of explicit chiral symmetry breaking, 
with an effective mass controlled by j. 

The solution of eq. (38) can reach a scale-invariant form. The only known fixed point 
of this equation corresponds to the Ising universality class, as only the fluctuations of the 
cr-field contribute to the evolution. The emergence of a critical theory around a = aj 
imposes two requirements: 

a) The first derivative of the potential around this point must vanish for A; — > 0. 

b) The fixed point of eq. (38) must be approached during the evolution. 

The solution of cqs.(38), (39) requires initial conditions. These are provided by the 
definition (33) for kx = T/6'2. (As before we use 62 = I for the calculations.) In particular 
we take 

Ckr = ^(0;T,^) = ^(a,;T,/i)-j (40) 
v',x - ^f^(..;T,;.)<^+i^(a,sT,;.)<5^ + l^(a,sT,;.)<5^ 



(41) 

and neglect the higher-order terms for small S. The first requirement for the emergence 
of a critical theory implies that c^j, must be chosen such that the solution of eq. (39) 
gives Co = 0. The second requirement can be satisfied if the initial condition (41) has two 
properties: The term ~ 5^ (the cubic term) vanishes, and the term ~ S (the mass term) 
has the appropriate critical value for the fixed point to be approached. The term ~ 5^ 
(the quartic coupling) need not be fine-tuned, as the Ising fixed point is characterized by 
one relevant parameter. For any given value of the quartic coupling, the fine-tuning of 
the mass term leads to the fixed point. For given j, we have three free parameters at 
our disposal: T, /i and the value around which we expand the potential. Therefore, 
we expect that the fixed point of eq. (38) can be approached for a unique choice of T, /i 
and aj. The values of these parameters mark the location of the critical endpoint on the 
phase diagram. 

In order to determine Cq, eq. (39) must be solved numerically. Alternatively, an 
analytical upper bound can be obtained through the relation 

In fiili-U -31.3 ^ 1"^-"^'+'' 



Ckr+J) aj Jkr 

The threshold function is monotonic and the maximum value of |Lf(u')| is |L'^(0)| = ^/^^. 
Thus, we obtain 



3 Tkr , 
S ■ (43) 



a 



For our choice of parameters (see below) the r.h.s. is significantly smaller than 1. This 
implies that there is only a small difference between cq and c^j,. As a result, the vanishing 
of Co implies \dUk^lda\ ((jj;T, /x) ~ j, according to eq. (40). 
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Figure 8: Determination of the initial conditions (40), (41) for the evolution equations 
(38), (39). 

The determination of the initial conditions for the solution of eqs. (38), (39) is shown 
in fig. 8. The first four derivatives of the potential Ukj. at the initial scale kx = T are 
plotted. The point (jj at which the third derivative vanishes is where we expect the critical 
theory to emerge. For the parameters of the model we have chosen (A = 0.1, /i^ = 2.3, 
K^/Pqr = 1.8), and i^/^/p^ = 0.72, T j ^fp^ ~ 0.127, we find aj^fp^ ~ 0.0858. We are 

considering an external source j / p^J^ ^ 1.12 x 10~^. The initial value of the parameter 
Ckj., defined in eq. (40), is Ckj./ p^^ ^ 1.74 x 10~^. We observe that [dUkj. / da]{aj]T, p) is 
close but not exactly equal to j. As a result [dVkj,/dS]{0]T, p) ^ 0. 

In fig. 9 we display elements of the solution of eqs. (38), (39). Wc observe that the 
parameter = [9Vfc/95](0; T, /x) starts from its initial value at the scale = T and 

1 /2 

effectively becomes zero at a scale k/ p^^ ~ exp(— 6). The rescaled mass of the pion fields 
ml{k)/k'^ = [{j + Cfe)/cTj]//c^ becomes much larger than 1 shghtly below this scale. As a 
result, the pions decouple completely from the evolution equation and Ck remains constant 
during the ensuing evolution. This guarantees that the first derivative of the potential at 
A; = vanishes at the point 5 = 0. 

The solution of eq. (38) approaches a fixed point in the infrared. This is demon- 
strated in fig. 9 by displaying: a) the minimum Sq — So/VkT of the rescaled po- 
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Figure 9: The evolution of: Ck = [dVk/d6]{0;T, fj,) , the rescaled a-field and pion masses, 
and the minimum of Vk{S) — {Vk{S;T, /i) — CkS) /{k^T). /i/y/poR — 0.72, T/ ■sJpqr ~ 

0.127, j/pl'^ ~ 1.12 X 10-^. 



tential Vk{5) = (Vfc(5;T,yu) — CkS) /{k^T); and b) the rescaled mass term w?„{k)/k'^ ~ 
[5^\4/9(5^]((5o; T, of the ci-field (or the shifted (5-field) at the minimum. Both these 
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quantities reach constant values at a scale k/ p^^ ^ exp(— 14). The subsequent evolution is 
analogous to the one studied in subsection 7.1. The potential, and the various parameters 
derived from it, stay constant near their fixed-point values for some "time" t — hi{k / p^^) ^ 
before finally deviating towards the symmetric or the broken phase. The renormalized 
theory at A; = exhibits a second-order phase transition. The parameter that has to be 
fine-tuned in order to go through the critical endpoint is a linear combination of T and 
p. 

The universahty class of the fixed point is determined by the number of fields with 
running masses comparable or smaller than k and the symmetries of the relevant part of 
the action. For k/ p\l^ ^ exp(— 6) only the cr-field fulfills this criterion. The pions have 
completely decoupled by the time the fixed point is approached. The scale-dependent 
potential v{5) has a Z2 symmetry that is preserved by the evolution. Crucial for this 
conclusion is the choice of a point Uj at which the cubic term of the potential U^^ vanishes. 
The elimination of the pion contribution from the evolution equation (38) through the 
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parametrization (37) is important as well. The pion fluctuations only affect the evolution 
of the term linear in 5 around aj (the parameter Cfe). This term evolves to zero and stays 
constant after the pion decoupling. These features can be achieved for specific values of 
T and /x that determine the location of the critical endpoint. As a result, the endpoint 
belongs to the Ising universality class. This can be confirmed by calculating the critical 
exponents, as was done in subsection 7.1. 

Before concluding this section, we should comment on the approximations employed 
in the study of the critical endpoint. The parametrization of eq. (37) is quite general. 
The derivation of the evolution equations (38), (39) made use of the crucial assumption 
v'i^{6) <^ Ck/vkFr. In the example we presented, the starting scale kx = T was sufficiently 
low for this relation to be satisfied to a good approximation during the whole evolution. 
We also neglected the terms higher than 6^ in the initial condition (41) for -y^^. Again, 
this is a reasonable approximation if we are interested only in a small region around aj. 

Going beyond these approximations in order to increase the precision of the study 
leads to significant technical complications. The reason is that the Z2 symmetry of Vk{S) 
is not guaranteed automatically if the approximations are relaxed. However, we believe 
that the essential conclusions remain valid. In order to have a critical theory at a point 
aj one must ensure that, ai k — 0, the first and third derivatives of the potential are zero 
there. These two conditions guarantee that aj would be a minimum, and that no second 
minimum would exist in its vicinity. The appropriate choice of aj and the fine-tuning of 
a linear combination of fi and T can achieve this goal. The fine-tuning of the orthogonal 
combination of n and T can be used in order to approach the fixed point that controls 
the critical theory. As the Ising fixed point is the only one known for the one-field theory 
after the pion decoupling, the resulting critical theory must be in the Ising universality 
class. 

The demonstration that this scenario is realized, if no approximations for the potential 
are made, must confront non-trivial technical complications. It is difficult to keep track 
simultaneously of the "non-universal" part of the potential corresponding to its first and 
third derivatives, and the "universal" part that approaches the fixed point. The reason 
is that the two parts scale differently with k. Our approximations in this work simplified 
this task by separating the evolution of from that of v'{S). The discussion of the general 
case merits a separate study. 

9 Summary and outlook 

The purpose of this study was to establish the appropriate framework for the analytical 
study of the QCD phase diagram, for the relatively low values of the chemical potential 
relevant for heavy-ion experiments. We employed an effective model of QCD, the linear 
quark-meson model. We considered only the case of two fiavours, as all the expected 
structure of the phase diagram appears already at this level. The model is similar to 
the (T-model of Gell-Mann and Levy [24] , with the two nucleons replaced by the u and d 
quarks. However, it can be extended in a straightforward manner in order to take into 
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account more flavours [25]. 

An efiicient framework must involve the renormalization group, in order to provide a 
reliable identification of the universality classes of the various second-order phase transi- 
tions. The Wilsonian (exact) formulation of the renormalization group permits also the 
detailed description of weak first-order phase transitions. We employed the formalism of 
the effective average action [16]. We concentrated on the evolution of the non-derivative 
part of the action, the mesonic potential Uk, as an effective infrared cutoff k is varied. 
The dependence of the potential on k is described by an evolution equation that takes the 
form of a partial differential equation. The potential at a scale kr equal to the tempeature 
T serves as the initial condition. For A; = the solution becomes equal to the effective 
potential, from which all physical information relevant for the vacuum structure can be 
extracted. 

We considered values of the parameters of the model that do not correspond to reahstic 
QCD. However, they were chosen such that the discussion could be carried out with the 
minimum use of numerics^. In particular, the potential at kx = T was derived analytically, 
using perturbation theory. Its form is determined by the potential at zero T and /i and 
the contributions from fiuctuations with characteristic momenta larger than /c^. The 
fermionic fiuctuations play a crucial role, as they can change the vacuum strucure of the 
potential by generating additional minima. They decouple at scales below T, because 
of the effective mass developed by the fermions at non-zero temperature. As a result, 
the entire fermionic contributions and the mesonic contributions with momenta larger 
than kr can be incorporated in the potential at /cy. The evolution below kx is effectively 
three-dimensional (dimensional reduction), as only the zero-modes of the mesonic fields 
are light at these scales. 

The other significant simplification afforded by our choice of couplings was the possibil- 
ity to neglect the wavefunction renormalization of both the fermionic and mesonic fields, 
and the evolution of the Yukawa coupling of the meson-fermion interaction. This meant 
that we could use a very simple truncation of the effective average action, which involves 
standard kinetic terms and a general scale-dependent potential for the a and pion fields of 
the mesonic sector. The evolution of the potential was determined by solving numerically 
a partial differential equation. As we mentioned above, the fermionic fluctuations affect 
only the initial condition for this equation. 

The phase diagram we derived has all the expected features. For zero current quark 
masses, there is a line of second-order phase transitions starting on the — axis. 
We confirmed that they belong to the 0(4) universality class, by calculating universal 
quantities such as the critical exponents /3 and v. For large there is a line of first-order 
phase transitions. 

The two lines meet at a tricritical point, with specific values (/x*, T^), where a second- 
order phase transition takes place, governed by the Gaussian fixed point. This results in 
mean-field behaviour, as we checked by calculating the relevant critical exponents. For 

"'^Thc parameters were also chosen such that the pattern of a second-order phase transition for fi = j = 
and increasing T and a first-order one for T = j = and increasing ij, appears. The reahstic values have 
been shown to predict this pattern [26, 27] . 
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II slightly smaller than we observed universal crossover behaviour (not to be confused 
with the analytical crossover) , as the initial influence of the Gaussian flxed point is slowly 
dominated by the more stable 0(4) fixed point near the critical temperature. For /i 
slightly larger than we observed very weak first-order phase transitions, for which the 
discontinuity in the order parameter approaches zero. 

We then took into account the effects of non-zero current quark masses by introducing 
a source term —ja. We discussed how the second-order phase transitions turn into an an- 
alj^ical crossover, while the first-order ones retain their qualitative character. We studied 
in detail the emergence of the critical endpoint of the line of first-order phase transitions. 
The pions are massive near the endpoint and we demonstrated their decoupling from the 
evolution equation. A second-order phase transition takes place at the endpoint, and we 
showed that it belongs to the Ising universality class. 

The study demonstrated that it is possible to obtain a complete and detailed picture 
of the QCD phase diagram using an effective theory such as the quark-meson model. The 
renormalization group is indispensable in the process. 

Realistic values for parameters of the model cannot be considered without having 
to rely on extensive numerical work. The most important issue concerns the value of 
the Yukawa coupling, which must significantly larger than 1 for the constituent quark 
masses to be reproduced correctly [29]. This implies that the evolution of this coupling 
and the wavefunction renormalization of the fields must be taken into account. Of par- 
ticular interest is the emergence of partial infrared fixed points at zero temperature and 
chemical potential [29, 27]. The derivation of the initial conditions for the effectively three- 
dimensional evolution at scales below T cannot be achieved analytically for large Yukawa 
couplings. However, it can be performed efficiently through the numericall integration of 
the appropriate evolution equations [27]. 

One technical issue related to the dimensional reduction of the theory at energy scales 
below T concerns the fermionic infrared cutoff that is used in the definition of the effective 
average action. A form must be devised that preserves the Lorentz structure of the kinetic 
term of free fermions for non-zero chemical potential, and in the same time guarantees 
fast decoupling of the temperature effects at low energy scales. 

Finally, the strange quark can be taken into account by considering the linear quark- 
meson model with an SU{3)]^ x SU{3)r chiral flavour group. The increase in complexity 
of the resulting evolution equations is signiflcant. However, the process can be carried 
out in a straightforward manner along the lines of refs. [25]. 

The present study has established that the linear quark-meson model is an effective 
theory of QCD that predicts all the expected qualitative structure of the phase diagram. 
The main aim of future studies along the lines highlighted above will be to determine 
with accuracy the location of the critical endpoint and the size of the area around it in 
which the critical behaviour persists. In this way, the possible effect of the endpoint on 
observable signatures in heavy-ion experiments will be determined quantitatively. 
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